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I.  Introduction 

This  Final  Report  is  a  summary  of  the  research  performed  under  the  grant 
N000 14-98- 1-02 12,  “Basic  Research  in  Thermoacoustic  Heat  Transport,”  which  was 
active  from  1  January  1998  to  31  May  2003.  The  grant  has  been  shared  between  the  two 
principal  investigators,  Robert  M.  Keolian  and  Anthony  A.  Atchley.  The  work  of  the  two 
groups  is  described  below  followed  by  a  combined  listing  of  external  output  in  the  form 
of  papers,  presentations  and  graduate  students  supported  at  least  partially  by  the  grant. 


Distribution  Unlimited 


II.  Keolian  Group 


A.  Objective 

The  guiding  objective  of  this  portion  of  the  grant  has  been  to  perform  basic 
research  that  lowers  the  barriers  to  Naval  or  commercial  adoption  of  thermoacoustic 
technology— namely  studying  issues  that  would  improve  the  energy  and  volume 
efficiency  of  engines  and  refrigerators. 


B.  Approach 

Our  technical  approach  has  been  in  several  directions:  1)  Model  a  thermoacoustic 
stack-based  waste  heat  driven  chiller  for  use  on  a  Navy  surface  combatant.  2)  Develop  a 
theory  for  thermoacoustic  refrigeration  and  sound  generation  in  standing  wave  devices 
that  do  not  use  a  stack.  We  call  these  “no-stack”  devices.  3)  Perform  experiments  on  a 
no-stack  device  which  models  a  radial  heat  driven  chiller  that  we  believe  is  compact 
enough  to  be  useful  to  the  Navy  and  commercially.  4)  Experimentally  map  out  the  heat 
transfer  and  drag  performance  of  various  heat  exchanger  types  in  the  presence  of 
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oscillating  flow.  5)  Construct  a  compact  thermoacoustic-Stirling  (traveling  wave  style) 
engine  for  use  in  electrical  power  generation. 


C.  Synopsis  of  Scientific  Accomplishments 

1.  Waste  Heat  Driven  Chiller 

Navy  destroyers  generate  substantial  amounts  of  waste  heat  in  the  form  of  the  hot 
exhaust  from  their  gas  turbine  engines  which  are  used  for  main  propulsion  and  electrical 
power  generation.  These  high  temperature  exhaust  streams  make  an  attractive  source  of 
high  quality  waste  heat  to  power  a  thermoacoustic  refrigeration  plant.  Such  a  plant  could 
provide  chilled  water  for  the  ship,  while  at  the  same  time  cool  the  exhaust  gases  and 
substantially  reduce  the  electrical  load  on  the  ship. 

A  thermoacoustic  refrigerator  uses  high  amplitude  sound  to  pump  heat  from  a 
cold  to  a  warmer  temperature.  A  thermoacoustic  engine  is  used  in  the  opposite  direction, 
to  generate  sound  from  heat.  What  is  described  here  is  a  heat  driven  thermoacoustic 
chiller  that  combines  both  types  of  devices.  Figure  1  shows  a  cut-away  sketch  of  a 
portion  of  a  waste  heat  driven  thermoacoustic  chiller  concept  that  we  have  modeled.  The 
complete  chiller  is  in  the  shape  of  an  annulus  that  surrounds  a  hot  exhaust  duct,  shown  by 
the  large  arrow.  Near  the  hot  inner  wall  of  the  annulus,  some  of  the  hot  exhaust  is 
directed  into  the  chiller  and  runs  through  a  hot  heat  exchanger  (labeled  Hot  HX),  drawn 
in  the  sketch  as  a  tube  and  fin  type,  although  other  types  of  exchangers  are  possible.  The 
next  item,  working  radially  outward,  is  a  stack,  shown  here  as  a  stack  of  washer  shaped 
rings.  On  the  outer  edge  of  this  inner  stack  is  another  heat  exchanger  kept  near  the 
ambient  temperature  by  an  ambient  temperature  fluid,  such  as  sea  water,  fresh  water 
cooled  by  sea  water,  or  outside  ambient  air.  Similarly,  towards  the  outside  of  the  annulus 
are  placed  a  cold  heat  exchanger,  an  outer  stack,  and  an  ambient  heat  exchanger.  Ducting 
above  and  below  the  heat  exchangers  is  provided  for  the  ambient  and  cold  heat  exchange 
fluid  streams.  In  the  annular  acoustic  space  surrounding  the  heat  exchangers  and  stacks 
is  the  working  gas,  which  can  be,  for  example,  a  pressurized  mixture  of  helium  and 
xenon,  for  the  best  performance,  or  air  at  ambient  pressure,  for  the  most  simplicity.  Pure 
helium,  or  helium-argon  mixtures  can  be  used  as  well  at  various  pressures,  all  giving 
various  tradeoffs  between  performance  and  simplicity. 

The  inner  stack,  with  its  hot  and  ambient  heat  exchangers,  is  the  source  of  sound. 
The  large  temperature  gradient  across  this  stack  causes  sound  to  spontaneously  start 
oscillating  in  the  radial  direction.  The  sound  shuttles  heat  out  of  the  hot  heat  exchanger, 
moves  it  down  the  temperature  gradient  along  the  inner  stack,  and  into  the  inner  ambient 
exchanger,  cooling  the  turbine  exhaust  gas  and  generating  sound  in  the  processes.  These 
radial  oscillations  of  the  working  gas  then  pump  heat  out  of  the  cold  heat  exchanger  on 
the  outer  side  of  the  annulus,  move  it  along  the  outer  stack  up  a  temperature  gradient,  and 
into  the  outer  ambient  heat  exchanger,  forming  a  refrigerator.  The  net  result  is  that  heat 
leaves  both  the  hot  exhaust  and  the  chilled  fluid  and  is  deposited  into  the  ambient  fluid, 
all  with  no  moving  parts  inside  the  annular  resonator. 
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Figure  1.  A  stack-based  waste  heat  driven  chiller  concept. 

Note  that  the  outer  surface  of  the  annular  chiller  is  kept  near  ambient  temperature 
by  the  outer  ambient  heat  exchanger.  One  can  think  of  the  chiller  as  being  like  an  “active 
lagging,”  adding  about  14  inches  to  the  radius  of  the  duct  in  one  design,  but  replacing 
part  of  the  lagging  normally  around  a  hot  duct,  thus  saving  space. 

Performance  of  the  chiller  is  dependent  on  the  type  of  gas  and  pressure  chosen  for 
the  working  fluid.  Computer  models  of  the  performance  were  made  for  an  annular  chiller 
surrounding  a  five  foot  diameter  exhaust  duct  carrying  combustion  products  at  980  °F, 
with  an  ambient  water  temperature  of  95  °F  and  a  chilled  water  temperature  of  45  °F,  for 
various  gas  and  pressure  combinations.  The  performance  ranged  from  8.8  tons  of  cooling 
(30.9  kW  removed  from  the  cold  heat  exchanger)  per  foot  of  active  duct  length  for  a 
helium-xenon  working  gas  mixture  at  3  atmospheres,  to  3.7  tons  of  cooling  per  foot  of 
active  duct  for  air  as  the  working  gas  at  1  atmosphere.  These  two  cases  pull  105  kW  and 
86  kW  of  energy  out  of  the  hot  exhaust  per  foot  of  active  duct,  respectively.  For  long 
active  ducts  the  performance  will  drop  as  the  temperature  of  the  exhaust  gas  is  reduced. 

Although  the  annular  geometry  of  this  chiller  achieves  good  volumetric 
efficiency,  the  energy  efficiency  of  stack-based  standing-wave  thermoacoustic  devices 
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such  as  this  is  respectable,  but  not  great— just  over  20%  of  the  ideal  Carnot  efficiency  at 
best.  It  was  thought  that  non-stellar  energy  efficiency  would  be  acceptable  in  a  waste 
heat  driven  system  where  the  input  energy  is  free.  It  should  be  noted  that  since  this  work 
was  performed  large  energy  efficiency  improvements  have  been  demonstrated  for 
thermoacoustic-Stirling  devices  by  our  colleagues  at  Los  Alamos,  as  described  in  item  5 
below. 


2.  No-stack  Theory 

In  a  conventional  standing  wave  thermoacoustic  device,  shown  in  Fig.  1  and  again 
in  Fig.  2(a),  about  half  of  the  total  dissipative  energy  loss  occurs  in  the  stack,  mostly  as  a 
result  of  the  irreversible  heat  transfer  between  the  working  gas  and  the  stack  with 
additional  losses  due  to  viscous  “scrubbing”  of  the  gas  moving  along  the  stack.  The 
function  of  the  stack  is  to  put  oscillating  fluid  elements,  shown  by  the  horizontal  arrows, 
thermally  in  series,  causing  thermal  energy  to  be  transferred  from  element  to  element. 
This  allows  a  thermoacoustic  refrigerator  or  engine  to  span  a  large  temperature  difference 
between  the  hot  and  cold  heat  exchangers  that  straddle  the  stack  while  using  a  modest 
acoustic  amplitude. 


(a) 


(b) 


|<—  Stack— >| 


Cold  Hot 

HX  HX 


Figure  2.  (a)  A  conventional  standing  wave  thermoacoustic  device  with  a  stack  sandwiched  between  hot 
and  cold  heat  exchangers  in  a  resonator,  (b)  A  no-stack  device  where  gas  elements  travel  across  the  no¬ 
stack  gap  directly  between  the  heat  exchangers. 


We  have  been  studying  an  alternative  system  [Wakeland  and  Keolian  (2002), 
“Thermoacoustics  with  Idealized  Heat  Exchangers  and  No  Stack”]  shown  in  Fig.  2(b), 
which  is  predicted  to  have  higher  efficiency  and  greater  simplicity  under  the  proper 
conditions.  The  stack  has  been  removed  and  the  acoustic  amplitude  is  much  increased  to 
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the  point  that  fluid  elements  are  able  to  travel  directly  between  the  hot  and  cold  heat 
exchangers.  It  is  the  “ram  jet”  of  thermoacoustic  devices. 

To  calculate  the  work  and  heat  flows  of  the  thermodynamic  cycle,  we  follow  the 
path  of  gas  elements  in  a  Lagrangian  way  between  the  two  exchangers.  Figure  3  shows 
the  case  of  heat  exchangers  which  are  assumed  to  be  perfect— a  gas  element  immediately 
comes  to  the  temperature  of  a  metal  fin,  shown  in  gray,  as  soon  as  it  comes  into  contact. 
Both  the  engine  cycle  (a)  and  refrigeration  cycle  (b)  are  shown.  Losses  are  calculated  in 
a  Eulerian  way,  accounting  for  thermal,  viscous  and  minor  (nonlinear)  losses  at  the  heat 
exchangers,  thermal-viscous  losses  at  the  resonator  walls,  and  thermal  conduction 
between  the  exchangers.  In  one  refrigerator  configuration,  with  the  assumption  of 
thermally  perfect  heat  exchangers,  the  model  predicted  that  the  efficiency  would  be  quite 
high,  at  38%  of  the  Carnot  efficiency. 


(a) 


(b) 


Figure  3.  Temperature  T  vs.  displacement  x  loops  for  no-stack  gas  elements  going  through  (a)  an  engine 
cycle  and  (b)  refrigeration  cycle.  Perfect  heat  exchange  is  assumed  between  the  gas  and  heat  exchangers, 
shown  gray,  at  a  hot  temperature  Th  and  cold  temperature  Tc.  Temperature  is  plotted  vertically,  Lagrangian 
position  of  the  element  is  plotted  horizontally. 


Figure  3  instead  shows  the  paths  of  gas  elements  for  the  case  of  imperfect  heat 
exchange.  The  open  area  of  the  loop  is  roughly  proportional  to  the  amount  of  heat  and 
work  transferred  by  the  element.  The  rate  of  heat  transfer  between  the  element  and  the 
heat  exchanger  fin  is  proportional  to  the  parameter  C=  1/(1000  rf),  where  ris  an  effective 
thermal  relaxation  time  between  the  gas  and  the  heat  exchanger  and/is  the  frequency. 
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x  [m]  *  [m] 


Figure  4.  The  temperature-position  paths  of  a  gas  element  in  the  case  of  imperfect  heat  exchange  for  a  no¬ 
stack  device.  The  rate  of  heat  exchange  between  gas  and  fin  is  proportional  to  the  parameter  C. 


We  found  [Wakeland  and  Keolian  (2003)  “Calculated  effects  of  pressure-driven 
temperature  oscillations  on  heat  exchangers  in  thermoacoustic  devices  with  and  without  a 
stack”]  that  pressure  swings  of  the  gas  within  a  heat  exchanger  cause  surprisingly  large 
average  heat  flows  even  in  heat  exchangers  that  are  nearly,  but  not  quite,  perfect.  So 
much  so  that  we  had  to  revise  the  efficiency  predictions  of  our  no-stack  paper  from  38% 
down  to  32%  of  the  Carnot  efficiency.  Thus  the  no-stack  approach,  when  optimized,  can 
be  expected  to  be  more  efficient  than  the  stack-based  approach,  but  not  as  good  as  the 
thermoacoustic- Stirling  devices  which  have  achieved  efficiencies  of  40%  of  Carnot. 
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3.  No-stack  experiments 

We  have  tested  some  of  our  no-stack  ideas  in  a  resonator  built  to  mimic  a  high 
power,  radial,  waste  heat  driven  chiller  similar  to  that  in  Fig.  1  but  without  a  stack.  We 
believe  this  geometry  is  compact  enough  for  Navy  ships  or  commercial  use.  We  built  a 
trapezoidal  resonator  to  model  a  small  section  of  an  annular  chiller.  The  resonator  holds 
a  no-stack  engine  at  its  narrow  end.  A  refrigeration  stage,  if  used,  would  go  at  the  wider 
end.  The  heat  for  the  engine  comes  from  electrical  heaters  that  substitute  for  the  hot 
combustion  products  of  a  gas  turbine. 


Figure  5.  No-stack  experiment.  The  experimental  apparatus  (on  the  right)  is  a  trapezoidal  resonator  with 
heat  input  (red  arrows)  on  the  narrow  end.  A  no-stack  engine  is  near  the  top  of  the  resonator.  The 
trapezoidal  geometry  is  motivated  by  being  a  small  section  of  an  annular  waste  heat  driven  chiller  concept 
(on  the  left). 


Much  to  our  surprise,  the  engine  achieved  onset,  but  ran  at  an  amplitude 
considerably  smaller  than  that  what  was  necessary  for  a  gas  parcel  to  travel  the  distance 
between  the  hot  and  cold  heat  exchangers.  We  have  made  a  numerical  model  of  these 
results  that  needed  to  include  the  heat  transfer  between  gas  parcels  in  the  acoustic 
direction  (an  effect  which  is  ignored  in  our  no- stack  paper).  These  experiments  are  not 
yet  complete. 
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4.  Oscillating  Flow  Heat  Transfer  Experiments 

By  placing  two  identical  heat  exchangers  in  an  oscillating  flow,  separated  by  a 
gap,  one  can  get  a  rather  straightforward  measurement  of  the  oscillatory  heat  transfer 
coefficient  under  various  conditions  by  measuring  the  heat  flow  between  the  exchangers 
and  dividing  by  half  of  the  temperature  difference  between  them.  The  drag  can  be 
determined  by  the  oscillatory  pressure  drop  of  the  flow.  These  measurements  are  useful 
as  input  to  our  no-stack  models,  but  should  also  be  useful  in  the  design  of  other 
thermoacoustic  devices  and  as  test  cases  for  validation  of  computational  acoustics 
calculations. 

Our  oscillating  flow  heat  exchange  and  drag  apparatus,  shown  in  Fig.  6,  relies  on 
similitude.  It's  experimentally  easier  to  build  scaled  up  models  of  heat  exchangers  and 
run  in  the  range  of  0.1  to  10  Hz  than  at  the  40  Hz  to  10  kHz  frequencies  normally 
associated  with  thermoacoustic  devices.  Larger  “acoustic”  displacements  and  penetration 
depths  allow  for  easier  and  more  accurate  construction  of  the  heat  exchanger  models. 
Pressure  and  thermal  instrumentation  yield  acoustic  loss  and  heat  transfer  information. 


Figure  6.  Oscillatory  flow  heat  exchange  and  drag  apparatus,  on  left,  and  experimental  heat  transfer  data 
(points)  and  correlations  (heavy  lines)  fitting  the  data,  on  right. 
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The  heat  transfer  measurements,  some  of  which  are  also  in  Fig.  6,  were  distilled 
into  an  effectiveness  model  for  heat  exchange  from  oscillatory  flow  between  parallel 
plate  heat  exchangers  [Wakeland  and  Keolian  (2003),  “Effectiveness  of  parallel-plate 
heat  exchangers  in  thermoacoustic  devices”].  The  drag  measurements  were  reduced  to 
fits  for  the  acoustic  resistance  of  parallel  plate  heat  exchangers  [Wakeland  and  Keolian 
(2003),  “Measurements  of  the  resistance  of  parallel-plate  heat  exchangers  to  oscillating 
flow  at  high  amplitudes”]  and  individual  screens  [Wakeland  and  Keolian  (2003), 
“Measurements  of  resistance  of  individual  square-mesh  screens  to  oscillating  flow  at  low 
and  intermediate  Reynolds  numbers”],  as  shown  in  Fig.  7.  We  believe  this  oscillating 
flow  data  and  their  fits  will  be  generally  useful  to  designers  of  thermoacoustic  devices, 
even  if  the  no-stack  concept  that  motivated  it  fades  from  use. 


Re 


hx 


Figure  7.  Oscillatory  flow  resistance  measurements  and  fits  for  parallel  plate  heat  exchangers,  on  left,  and 
individual  wire  screens,  on  right. 


5.  Compact  Thermoacoustic-Stirling  Engine 

Scott  Backhaus  and  Greg  Swift  of  the  Los  Alamos  National  Laboratory  have 
vastly  improved  the  possible  efficiency  of  thermoacoustic  devices  with  their  invention  of 
the  thermoacoustic  Stirling  prime  mover  (Backhaus  and  Swift,  J.  Acoust.  Soc.  Am.  107, 
3148-3166  (2000)).  However,  their  engine  is  quite  large  and  has  an  unwieldy  shape 
which  seems  unsuitable  for  placement  on  a  Navy  ship.  We  have  constructed  and  are 
testing  a  more  compact  version  of  their  engine  that  will  be  used  in  future  work  for  the 
purpose  of  generating  electrical  power.  We  have  given  up  some  of  the  niceties  of  the 
flow  present  in  the  Backhaus  and  Swift  apparatus  in  exchange  for  compactness,  hoping 
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that  the  efficiency  hit  we  expect  to  take  will  be  more  than  made  up  by  the  inherent 
superiority  of  the  Stirling  cycle  compared  to  the  conventional  stack-based  thermoacoustic 
cycle.  These  experiments  are  ongoing. 


Figure  8.  Compact  thermoacoustic-Stirling  heat  engine. 


III.  Atchley  Group 


A.  Objective 

The  overall  technical  objective  of  this  portion  of  the  grant  is  to  conduct  a  careful 
experimental  investigation  of  high  amplitude  acoustic  phenomena  in  geometrically- 
simple  resonators.  Applications  of  high  amplitude  acoustics,  such  as  thermoacoustics, 
have  stimulated  the  search  for  a  fuller  understanding  of  high-amplitude  acoustics  in 
closed  systems  than  is  accessible  through  analytical  theory.  As  a  result,  several 
computational  approaches  have  been  investigated.  While  the  ultimate  goal  of  such 
modeling  is  application  to  practical  configurations,  the  first  steps  are  validation  against 
experimental  results  acquired  in  geometrically-simple  systems.  The  purpose  of  this 
portion  of  our  research  is  to  aid  the  development  of  these  models  by  providing 
measurements  of  various  nonlinear  acoustic  phenomena  in  computational  tractable 
resonator  geometries  against  which  to  compare  numerical  predications. 
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B.  Approach 

The  technical  approach  has  been  to  focus  on  two  aspects  of  the  problem:  1) 
Rayleigh  streaming  in  cylindrical  resonators,  2)  the  time-averaged  pressure  change  across 
an  abrupt  change  in  resonator  cross  section. 


C.  Synopsis  of  Scientific  Accomplishments 

1.  Rayleigh  Streaming  in  Cylindrical  Resonators 

Classical  theories  of  Rayleigh  streaming  in  cylindrical  tubes  predicts  the  flow  in 
the  core  outside  the  boundary  layer  to  be  a  quadratic  function  of  radial  position. 
However,  these  theories  are  applicable  to  very  low  Reynolds  numbers,  consistent  with 
axial  components  of  the  streaming  velocity  of  a  few  mm/s  or  less.  At  higher  Reynolds 
numbers,  inertial  effects,  which  are  neglected  in  classical  theories,  become  important.  As 
pointed  out  by  Menguy  and  Gilbert  (. Acustica ,  86  249-259,  2000),  one  consequence  of 
including  inertial  effects  is  that  the  profiles  are  no  longer  parabolic.  We  have  used  laser 
Doppler  anemometry  to  measure  the  time-averaged  flow  in  a  cylindrical  resonator.  Our 
measurements  show  a  non-parabolic  profile,  as  can  be  seen  in  Fig.  8.  However,  the 
deviation  from  parabolic  is  much  greater  than  can  be  explained  by  including  inertial 
effects.  Measurements  also  show  that  the  streaming  velocity  starts  off  near  the  classical- 
theory  value  and  evolves  over  a  time  scale  of  several  minutes.  This  evolution,  which  also 
has  never  been  reported  before,  is  correlated  with  the  small  acoustically  generated 
temperature  gradients  induced  along  the  resonator  wall.  These  results  imply  that 
computational  models  need  to  account  for  long-term  dynamics  and  thermoacoustic  heat 
exchange  with  the  wall.  Details  of  these  measurements  are  provided  in  the  manuscript 
M.  W.  Thompson  and  A.  A.  Atchley,  “Optical  measurement  of  Rayleigh  streaming  in  a 
standing  wave  with  a  temperature  gradient,”  submitted  to  J.  Acoust.  Soc.  Am.,  June  2003. 
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radial  position,  r 


Figure  9.  Axial  streaming  velocity  across  a  cylindrical  tube.  “Controlled  gradient”  data  has  water  circulating 
around  the  outside  of  the  tube,  keeping  the  tube  at  nearly  constant  temperature.  “Uncontrolled  gradient” 
data  is  taken  without  the  circulating  water. 


2.  Time  Averaged  Pressure  Drops  at  Transitions 

Abrupt  changes  in  cross  section  occur  in  thermoacoustic  engines  at  junctions 
between  the  resonator  and  heat  exchangers,  between  heat  exchangers  and  the  stack,  at  jet 
pumps,  and  in  the  resonator  itself.  These  changes  in  cross  section  can  have  both 
dissipative  (resistive)  and  non-dissipative  (reactive)  effects,  influencing  both  the  nuisance 
loss  in  the  resonator  and  its  operating  frequency.  Unnecessary  dissipation  leads  to  a 
reduction  in  efficiency  of  an  acoustic  engine.  Unaccounted-for  reactance  can  lead  to  an 
engine  operating  at  a  non-optimal  frequency,  also  reducing  its  efficiency.  The  influence 
of  a  change  in  cross  section  depends  upon  the  geometrical  properties  of  the  change,  its 
location  in  the  sound  field,  and  the  amplitude  of  the  sound  field.  As  is  the  case  with 
streaming,  a  thorough  understanding  of  both  the  good  and  bad  consequences  of  acoustic 
flows  through  changes  in  cross  section  is  important  for  designing  acoustics  engines. 

We  have  developed  an  experimental  apparatus,  shown  in  Fig.  9,  that  allows  for 
temperature-controlled  measurement  of  the  acoustic  power  delivered  to  a  resonator  as 
well  as  the  absolute,  time-dependent  pressure  at  several  locations  along  the  resonator 
wall.  Measurements  of  the  time-averaged  pressure  at  the  narrow  end  of  a  stepped- 
resonator  are  in  very  good  agreement  with  computational  results  by  Sparrow,  et.  al.,  (also 
ONR  sponsored),  based  on  a  two  dimensional  unsteady  compressible  high-order  Navier- 
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Stokes  code  specialized  for  aeroacoustics  problems.  The  agreement  demonstrates  that  the 
types  of  measurements  pursued  here  are  capable  of  serving  their  intended  function,  to 
validate  computational  models.  This  part  of  the  project  has  been  the  subject  of  several 
presentations  at  meetings  of  the  Acoustical  Society  of  America.  The  details  will  be  found 
in  A.  Doller’s  dissertation,  scheduled  to  be  completed  at  the  end  of  December  2003,  and 
subsequent  publications. 


Figure  10.  Apparatus  for  measurement  of  mean  pressure  differences,  with  comparison  between 
experimental  data,  computational  results,  and  Bernoulli  theory. 
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